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Abstract 

We give a new computational method to obtain symmetries of ordinary 
differential equations. The proposed approach appears as an extension of 
a recent algorithm to compute variational symmetries of optimal control 
problems [Comput. Methods Appl. Math. 5 (2005), no. 4, pp. 387-409], 
and is based on the resolution of a first order linear PDE that arises as a 
necessary and sufficient condition of invariance for abnormal optimal con- 
trol problems. A computer algebra procedure is developed, which permits 
to obtain ODE symmetries by the proposed method. Examples are given, 
and results compared with those obtained by previous available methods. 
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1 Introduction 

Soplius Lie was the first to introduce the use of symmetries into the study of 
differential equations, Emmy Noether the first to recognize the important role 
of symmetries in the calculus of variations. Currently, all the computer algebra 
systems that address differential equations provide several tools to help the user 
with the analysis of Lie symmetries. Recently, the authors developed a computer 
algebra package for the automatic computation of Noether variational symme- 
tries in the calculus of variations and optimal control [5J , now available as part of 

the Maple Application Center at http : //www.maplesof t . com/applications/app_center_view. aspx?AID=19 

The omnipresent tools for Lie symmetries provide a great help for the search 
of solutions of ODEs, their classification, order reduction, proof of integrability, 
or in the construction of first integrals. From the mathematical point of view, a 
ODE symmetry is described by a group of transformations that keep the ordi- 
nary differential equation invariant. Depending on the type of transformations 
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one is considering, different symmetries are obtained. An important class of 
symmetries is obtained considering a one-parameter family of transformations, 
which form a local Lie group. Those transformations are often represented by 
a set of functions known as the infinitesimal generators. From the practical 
point of view, the determination of the infinitesimal generators that define a 
symmetry for a given ODE is, in general, a complex task [HI UT] . To address the 
problem, we follow a different approach. 

We propose a new method for computing symmetries of ODEs by using 
a Noetherian perspective. Making use of our previous algorithm [S], that has 
shown up good results for the computation of Noether variational symmetries of 
problems of the calculus of variations and optimal control, we look to an ODE 
as being the control system of an optimal control problem. Then, we obtain 
symmetries for the ODE by computing the abnormal variational symmetries of 
the associated optimal control problem. 

This paper is organized as follows. In ij2l the necessary concepts associated 
with variational symmetries in optimal control are reviewed. The new method 
for computing symmetries of ODEs is explained in The method is illustrated 
in [21 where we compute symmetries for three distinct ODEs and compare the 
results with the ones obtained by the standard procedures available in Maple. 
We end the paper with some conclusions and comments SJSl The definitions of 
the new Maple procedure that implements our method are given in Appendix. 



2 Symmetries in optimal control 

Without loss of generality, we consider the optimal control problem in Lagrange 
form: to minimize an integral functional 

/[x(.),u(.)] - f L{t,x{t),u{t))dt (1) 

J a 

subject to a control system described by a system of ordinary differential equa- 
tions of the form 

±{t) ^ ip{t,Mt),n{t)) , (2) 

together with appropriate boundary conditions. The Lagrangian L : M. x R" x 
^ M and the velocity vector : R x R" x R™ ^ R" are assumed to be 
continuously differentiable functions with respect to all their arguments. The 
controls u : [a,b] R™ are piecewise continuous functions; the state variables 
x: [a,b] R" continuously differentiable functions. 

The celebrated Pontryagin Maximum Principle [TU] (PMP for short) gives 
a first-order necessary optimality condition. The PMP can be proved from a 
general Lagrange multiplier theorem. One introduces the Hamiltonian function 

H{t, X, u, ipo, -tp) = il^QL{t, X, u) -f ■0'^ ■ ip{t, X, u) , (3) 

where (i/'o, ''/'(•)) ^i'^ the "Lagrange multipHers", with ?Ao < a constant and 
a rt-vectorial piecewise C-'^-smooth function, and the multiplier theorem 
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asserts that the optimal control problem is equivalent to the maximization of 
the augmented functional 



J[x(.), u(.), ^0, ^(O] = / {H{t, x(t), u{t), tlJo,i>{t)) - ^{t)'^ ■ x(i)) dt . (4) 

Definition 1. A quadruple (x(-), u(-), i/jqi ''/'(•)) satisfying the Pontryagin Max- 
imum Principle is said to be a (Pontryagin) extremal. An extremal is said to 
be normal when ipo ^ 0, abnormal when tpo — 0. 



Let : [a, 5] X M" X R" X R X M" ^ R X : 
group of transformations of the form 



be a one-parameter 



h''(t,x, u, ■00,1/') = 

{h\{t, X, u, Vo, lA), hj(i, X, u, V'o, i^), hu(i, x, u, -00, ■0), h^(i, x, u, V'o, ■0)) • (5) 

Without loss of generality, we assume that the identity transformation of the 
group ([5]) is obtained when the parameter s is zero: 

/ij (t, X, u, -00, lA) — t: hx(i, X, u, 007 ■0) = X, 
h°(t,x,u,0o,'0) u, h^(t,x,u,0o,'0) = ip. 

Associated with a one-parameter group of transformations ([5]) , we introduce its 
infinitesimal generators: 



d 

T(i,x,u,0o,i/') = 



U(t,x,u,'0o,i/') TphJ 

OS 



s=0 



d 

, X(t,x,u,0o,'0) = — 
d 

, *(t,x,u,0o,'0) - ^h: 



(6) 



s=0 



We can define variational invariance using the augmented functional ^ and 
the one-parameter group of transformations ([5]) or an equivalent condition in 
terms of the generators (|6]): 

Definition 2 ([3lll3j). An optimal control problem (Ip-fl|j is said to be invariant 
under ([6]) or, equivalently, ([6]) is said to be a symmetry of the problem (Op-iOi 

V 



at ox 



U+ — ■* 

ovl dtp 



dX 



AT 



with H the Hamiltonian 



A computational algorithm to obtain the infinitesimal generators T, X, U, 
and 'S' that form a variational symmetry ([7]) for a given optimal control problem 
([iJ-Q was developed in [5]. Here we remark that the abnormal variational 
symmetries (i.e. the ones associated with tAq — 0) obtained by the method 
introduced in [5] provide symmetries for ordinary differential equations. 
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3 Computing ODE symmetries from an optimal 
control perspective 

To ODE symmetries (see e.g. [3 [TS]) we associate a smaller group of infinites- 
imal transformations than that used in optimal control: we are only interested 
in transformations of the independent variable t and transformations of the de- 
pendent variables x. Let : R x R" ^ M x R" be a one-parameter group of 
transformations of the form 

g^(t,x) = (5,^(t,x),gj(i,x)), 

with g*'(i,x) = (i,x;). Let us denote the respective infinitesimal generators by 



d 



(8) 



Our method begins by identifying a richer set of variational symmetries in the 
form ([S]), from which we then obtain ([5]). In this section, we explain in detail 
how to arrive to the set of generators ([5]) that keep an ODE invariant. 

We are interested in determining symmetries for systems of ODEs in the 
canonical form 



(''") A /'+ • ('•1-1) 



1 Vn J Vn ■ ■ ■ Vn ) 



1 ?/n J 2/ti • ■ • Vn j 



(9) 



where functions 0^ : R x R^"=i''' R, /c = 1, . . . , n, are continuously differen- 
tiable with respect to all their arguments. To write the system ^ of differential 
equations as a control system ^ , we begin by converting it as a system of equa- 
tions of first order. For that we introduce a new set of variables, represented by 
the vector x: 



X = [Xi,... 

= yi,yi, 

where r = T,"^^ri. With this notation, we get the control system 



jVi I ■ ■ ■ I ym ym ■ ■ ■ J Vn 



(10) 



Xi 



Xri — 1 



X2 



= 01 , 



, Xri + ...+r„ 



tTi+...+r„ I 



0n(i,x) , 



(11) 



with r state variables but no control variables (i.e., (fTTj) is a particular case of 
where (p does not depend on u). 
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To use the formalism of optimal control and the notion of variational sym- 
metry [5] , one thing is missing: the existence of an integral functional H]) to be 
minimized, and whose Lagrangian L enters into the definition of the Hamilto- 
nian thus being necessary for computing symmetries by ([7]). However, if we 
restrict ourselves to the abnormal case, where the cost functional has no role, i.e. 
if we fix t/iq = 0, then the Hamiltonian H does not depend on the Lagrangian L 
and we can look to our system (jlip as an optimal control problem. If we only 
consider the abnormal case, the control system ([TT|) is everything one needs to 
write ([7]) and find symmetries. 

We are now in conditions to use our Maple optimal control package [5J and 
its procedure Symmetry to obtain symmetries for systems of ODEs We 
only need to rewrite ([9]) as in (fTTj) and then call function Symmetry of [5] for the 
abnormal case. After using this technique with several concrete examples, and to 
be able to compare the obtained results with the ones from standard techniques, 
we concluded that a great manual effort is necessary at each particular problem 
in converting the initial system into the canonical form, then to (jlip . and finally 
recovering the initial notation to compare the results with those obtained by the 
tools already available in the Computer Algebra System Maple. To do all the 
process in an entirely automatic way, and also to optimize the algorithm, we 
define here a new Maple function odeSymm (see Appendix) whose purpose is to 
compute symmetries for systems © of ODEs. 

The algorithm 

We consider the abnormal Hamiltonian 

i/(t,x,V') =^^-'^(t,x), 
where the velocity vector is given by 

if{t,yi) = [a;2,...,a;ri,(^i(t,x),Xri+2,...,a;ri+r2>2(i,x), ••• 

• • ■ , Xj.^_|_____|_^^_^_|_2, . . . , ^r; ^n(^j x)] 

Condition ([7]) simplifies to 

Thus, given a system of ODEs, we determine the infinitesimal generators ^ and 
T] ([5]), which define a symmetry of the given ODEs, in the following way: 

1. First we rewrite the given system of ODEs in the form ([9])0 

2. We represent the dependent variables yi and their derivatives i = 
1, . . . , n, J = 1, . . . , 7'i — 1, present in functions 0^ of the canonical system 
([9]), by a new set of variables x, in accordance with pO|) : 

^We only deal with differential equations that can be written in the canonical form l(9] l. 
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(12) 
(13) 



3. We define the velocity vector ([T^ : 

4. We substitute ip and its partial derivatives into 

5. From equation (jl3p . we determine the variational generators T(t,:x.,xp), 
X(t, X, V') and \I>(t,x, t/j); 

6. In the results obtained, we go back to the initial notation of variables, by 
means of the inverse relations of (1101): 



7. From the obtained set of variational generators T, X and '3/, we extract 
the subset of generators ^ and rj, 

^ = T, r]i=X^y^^-i z = l,...,n, 

' fc= 1 

which represent symmetries for the given system of ODEs. 

The infinitesimal generators obtained in step [5] are functions of the auxiliary 
variables x. Since in step [6] the variables x resume to its initial meaning, we 
conclude that our method is able to give dynamic symmetries (cf. Example [21). 
Indeed, the generators may involve derivatives of the dependent variables: 

(e,r/) = (e(t,y,y,...),r7(i,y,y,---)) 

with (<,y,y, ...) = {t,yi,yi . . .y['''~^\ ■ ■ ■ , y„, y„ ■ • ■ 

We now address the non-trivial part of our seven-step algorithm, which re- 
sides precisely in step [5j the determination of the associated variational gener- 
ators. We use the following strategy. Expanding the total derivatives 

x uip 
we write equation (|13p as a polynomial 

A{t, X, xp) + B{t, X, • X -I- C(t, X, • i/> = (14) 
in the 2r derivatives x and xp: 

,T [dip dip dT ax\ 

+ — -t/^^-— -x 15 

ox ox / 

The terms in (1151) . which involve derivatives with respect to vectors, are ex- 
panded in row-vectors or in matrices, depending, respectively, if the function is 



dT 


dT 


dT . 


dT 


dt 




dx 


dip 


dX 


dX 


dX . 


dX 


"d^ 


dt " 


dx 


^ dip 
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a scalar function or a vectorial one. Equation fTS]) is a differential equation in 
the 2r + 1 unknown functions T, Xi, . . . , Xr and tf'i, . . . , ^Z'r- This equation 
must hold for all ii, . . . , Xr, ^i, . . . , ipr, and therefore the coefficients A, B, 
and C of polynomial (I14p must vanish, i.e, 



_*T^^T.^.^_^T.^^O^ (16) 

ax ax 
,T 9r ,T 9X 
O'lp dtp 

Although a system of partial differential equations, solving (|16p is possible using 
the Maple command pdsolve because the system is of the first order, homoge- 
neous, and linear with respect to the unknown functions and their derivatives. 
We also remark that since system is homogeneous, we always have, as 
trivial solution, (T, X, *) = 0. 

When dealing with ODEs with several dependent variables and high-order 
derivatives, the number of calculations to be done is big enough, and the help 
of the computer is more than welcome. We use the computer algebra system 
Maple 10 to define a new procedure odeSymm that does all the cumbersome 
computations for us - all the steps [T] to [7] of our algorithm. 

Our procedure odeSymm receives, as input, a system of ODEs, and returns, 
as output, a family of symmetries (^,J7) - see definition of procedure odeSymm 
in Appendix. To optimize the resolution of p6p . we give the possibility to 
pass several optional parameters to odeSymm. These optional parameters are 
described in the Appendix and illustrated with concrete examples in fjH Here 
we just mention that, by default, we use the method of separation of variables 
(see [SI [15) to solve (fTG]) . More precisely, we follow (21: the generators are 
replaced by the sum of unknown functions, one for each variable. For example, 
T{t, xi,X2) = Ti{t) -f T2{xi) + T^{x2)- Through the optional parameters, one 
can use the default solving process of the Maple solver pdsolve or other specific 
methods (cf. Example [T]). 



4 Illustrative examples 

To show the functionality and the usefulness of our new procedure odeSymm, we 
consider three concrete problems found in the literature. All the examples were 
carried out with Maple version 10 on a 1.4GHz 512MB RAM Pentium Cen- 
trino. The running time of procedure odeSymm is indicated, for each example, 
in seconds. 

Example 1 (Kamke's ODE 120). We begin with a first order ODE found in 
Kamke's book 

> ode:= t*diff (y(t) ,t)-y(t)*(t*ln(t"2/y(t))+2)=0; 
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orfe:=t|,(t)-,W (tin (1^) +2)^0 

To obtain symmetries of the equation we use our Maple procedure odeSymm with 
the additional parameter hint=noint. This means that we will use the default 
method of resolution of PDEs of the Maple solver pdsolve. If the optional 
parameter hint is not used (see Examples\^and\^below) , our procedure odeSymm 
uses the method of separation of variables. We obtain the following infinitesimal 
generators (0.72 sec): 

> gerad:= odeSyimn(ode , y(t) , split, hint=nohint) ; 



gerad :— 



, 1 y 



One can test the validity of the obtained symmetries with the symtest command 
of the DEtools Maple package: 

> map (DEtools [symtest] , [gerad], ode, y(t)); 

[0,0] 

The symtest confirm that the infinitesimal generators leave the given ODE in- 
variant, i.e., the generators obtained by our method give indeed a symmetry to 
Kamke's ODE 120. It is interesting to remark that, without the knowledge of 
the computed symmetries, the ODE Maple solver dsolve is not able to integrate 
the ODE: 

> dsolveCode, y(t) , class); 

However, when one gives to the Maple solver the infinitesimal generators found 
by our method, the ODE is correctly solved: 

> dsolveCode, y(t) , HINT= [gerad] ) ; 

It is also interesting to note that our method is able to find one symmetry that 
is different from the ones obtained using the standard methods of the literature. 
The Map I e system provides nine different algorithms to compute symmetries of 
ODEs through the command symgen of the DEtools package. All the available 
schemes for determining the infinitesimal generators - option way = all - are 
not able to identify our pair of infinitesimals [C = Oj ^ = •' 

> DEtools [symgen] (ode , y(t) , way=all) ; 



^ = 0,r? = yln — 

y 



Example 2 (Damped Harmonic Oscillator). We consider a harmonic oscillator 
with restoring force —kx, emerged in a liquid in such a way that the motion of 
the mass m is damped by a force proportional to its velocity. Using Newton's 
second law one obtains, as the equation of motion, the following second order 
differential equation |2l PP. 432-434]-' 

> EL:= m*diff (x(t) ,t,t)+a*diff (x(t) ,t)+k*x(t)=0; 
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EL ■- m—xit) + a—x(t) + kxit) = 
dt^ at 

The symmetries for this equation are easily obtained with our Map I e procedure 
odeSymm (1.21 sec) 

> gerad:= odeSyimn(EL, x(t), split); 



mx 



at t^Jo?^ihr, 



_ at t^a^-ikr, 



One can confirm that these infinitesimals represent valid symmetries for the 
differential equation: 

> map(DEtools [symtest] , [gerad] , EL, x(t)); 

[0,0,0,0,0] 

Note that the output of our odeSymm procedure includes a dynamical symmetry: 
the derivative of the dependent variable is present in the second pair of obtained 
infinitesimal generators. 

Example 3 (Kepler's problem). We now consider the Kepler's problem: a prob- 
lem of the calculus of variations - see 114\ P- 21 7]. In this case, the Lagrangian 
depends on two dependent variables qi and 52 ■' 

T / ■ \ rn , 2 2\ K 

Lit, q, q) = ^ (qi +q2) + , „ „ ■ 
^ V<ii + % 

We will use the proposed method to determine symmetries for the corresponding 
Euler- Lagrange differential equation. The Euler- Lagrange equation is trivially 
obtained using our package of the calculus of variations Example 5.2]: 

> L:= m/2*(v[l] ~2+v[2] ~2)+K/sqrt(q[l] ~2+q[2] -2) ; 

L ■= - m{vi +V2 ) + 



> EL:= CLawslCV] [EulerLagrange] (L, t, [q[l],q[2]], [v[l] ,v[2]] ) ; 



-m — q2{t) — = } 

In this case, the Euler- Lagrange equation is a system of two second order ODEs. 
Our odeSymm procedure is able to determine symmetries for systems of differen- 
tial equations as well ('13.32 sec): 

> odeSyiran(EL, [q[l] (t) ,q[2] (t)] , split); 
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[e = 0, r?i 



-12, m 



-t, TJl = qi, TJ2 = (72 



,[^ = 1,7^1= 0, r?2 = 0] 



It is worth to mention that this example can not be handled by the algorithms 
available in Maple. Indeed, the Maple command symgen that looks for a symme- 
try generator for a given ODE is not able to address more than one dependent 
variable. 



5 Conclusions 

We have used the CAS Maple to define a new computational procedure that de- 
termines, in an automatic way, symmetries of ODEs. The automatic calculation 
of symmetries is a subject much studied under the theory of differential equa- 
tions, with many results and applications in many different areas. Our main 
novelty is the presentation of a new algorithm, alternative to existing ones, which 
looks to symmetries of ODEs as particular cases of Noether- variational symme- 
tries. As explained in iJ31 our algorithm involves the resolution of a first order, 
homogeneous, and linear PDE, which is the abnormal case of the necessary and 
sufficient condition of invariance for problems of optimal control studied with 
Noether's theorem [5l[12]. Interesting points of the proposed method are: (i) it 
is based on a new approach to the subject - in particular, it is different from all 
the nine alternative algorithms available in Maple; (ii) allows us to obtain dy- 
namic symmetries for ODEs of any order; (iii) allows to determine symmetries 
for systems of ODEs, when the analog simgen Maple command of the DEtools 
package can only obtain solutions for a single ODE. 
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Appendix: the new Maple procedure odeSymm 

The procedure odeSymm, introduced in this paper, has been implemented for the 
computer algebra system Maple (version 10). The complete Maple definitions 



can be freely obtained from |http : / /www . ipb . pt/ "^gouveia/ odeSynun . htm to- 



gether with an online help database for the Maple system. 
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odeSymm computes the infinitesimal generators which define the symmetries 
of the ODE, or system of ODEs, specified in the input. As explained 
in section [31 this procedure involves the resolution of a system of partial 
differential equations. We have used the Maple solver pdsolve, using, as 
preferential method, the separation of the variables by sum. 

Output: 

- one or more lists of symmetry generators for a given ODE, or system 
of ODEs ([^ =?, 7?i =?, rj2 =-!,..., Vn =%■■■). 

Syntax: 

- odeSymm(ode, x(t), opts) 

Input: 

ode - ordinary differential equation, or a set or list of ODEs; 

x(t) - any indeterminate function of one variable, or a list of them, repre- 
senting the unknowns of the ODE problem; 

opts - (optional) specify options for the odeSymm command, where opts is 
one or more of the following: 

allconst - When this argument is given, the output presents all 
the constants given by the Maple command pdsolve. By default, 
that is, without option allconst, we eliminate redundant con- 
stants; this is done by our Maple procedure reduzConst, which 
is a technical routine, and thus not provided here. Essentially, 
the procedure transforms in one constant each sum of constants. 
The interested reader can find the Maple file with its definition 
at http : //www . ipb . pt/~pgouveia/odeSymm . htm, 

mindep - When one wants to restrict to the minimum the dependen- 
cies of the infinitesimal generators: ^{t) and rj(x). By default, 
that is, in the absence of options mindep and alldep, the follow- 
ing dependencies are considered: S,(t) and ?7(t,x); 

alldep - All possible dependencies for the infinitesimal generators: 
£_{t, x, i/)) and T]{t, x, ip); 

split - When this argument is given, the procedure invoke the split 
command to divide the resultant set of infinitesimal generators 
into uncoupled subsets, by fixing the values for all the constants 
given by the Maple command pdsolve. The procedure split is 
a technical routine. The interested reader can find the Maple file 
with its definition at http : //www . ipb . pt/'^gouveia/ odeSymm. htm^, 

showdep - Shows, in the obtained solution, all the dependencies of 
the generators; otherwise only the name of the generators is 
shown; 
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showt - Shows, in the obtained solution, the dependence on the time 

variable (independent variable); otherwise, the time variable is 

omitted as a function parameter; 
showgen - Shows, in the obtained solution, besides the infinitesimal 

generators ^ and rj, the augmented set of variational generators, 

T, X and *; 

hint=<i/alMe> - Indicate a method of solution of the PDE system 
(jl6p . where <value> is one of "+V*\ or any other expression 
allowed by the command pdsolve, being also possible to use 
hint=nohint for the case one wants to use the standard method 
of resolution of Maple; by default, the system is solved by sepa- 
rating the variables by sum (hint= 

odeSymm : = proc (ODEs : K ' = ' , list (' = '), set (' = ')}■, 

depvars : : {function, list (function) }) 
local n, tt , XX, pp, k, vX, vPSI, syseqd, sol, IstGerad, valGerad, phi, 
vphi, Ipsi, vpsi,Hi, t, Sr, xO, r, aux, mapx, sys , xieta, sol2; 
unprotect (Psi) ; unassign( 'T' ) ; unassign( 'X' ) ; unassign( 'Psi ' ) ; 
unassignC 'psi ' ) ; 

Hi : =subs (select (type , [args [3 . . -1] ] , '= ' ) ,hint) ; 
if Hi='hint' then Hi:='+'; fi; 
n:=nops (depvars) ; 

if n=l then xO : = [depvars] else xO:=depvars fi; 
t:=op(l,xO[l]) ; 
r: = [] : 

for aux in xO do 

for k from 1 by 1 while evalb(subs (diff (aux, t$k)=_zzz , DDEs) OODEs) do 
od; 

r: = [r[] , k-1] : 
od: 

Sr : =sum(r [' i '] , ' i ' =l..n); 

mapx : = [seq(xO [i] =_x [l+(sum(r ['k'] , 'k' =1 ..i-1))], i=l..n)]; 

mapx : = [mapx [] , seq(seq(dif f (xO [i] , t$j ) =_x [j + l+suiii(r [ 'k' ] , 'k'=l. .i-1)] , 

j=l. .r[i]-l),i=l. .n)] ; 
mapx : = [mapx [] , seq(dif f (xO [i] ,t$r [i] )=_xx [i] ,i=l. .n)] ; 
mapx :=[seq(mapx[nops (mapx) +l-i] ,i=l. . nops (mapx) )] ; 
sys : =subs (mapx, DDEs) ; 
if n=l then solve (sys , {_xx [1] }■) 

else solve ({sys []}, {seq(_xx [i] , i=l . .n) }) fi; 
phi : =subs (7. , [seq(_xx [i] , i=l . .n)] ) ; 

vphi : =Vector ( [seq( [seq(_x [j] , j=2+sum(r ['k' ] ,'k'=l..i-l).. sum(r [ 'k'] , 

'k'=l..i)), phi[i]][],i=l..n)]); 

xO:= [seq(_x[i], i = 1 .. Sr)] ; 

if Sr>l then Ipsi := [seq (psi [i] , 1=1 .. Sr)] else lpsi:=[psi] fi: 

vpsi :=Vector [row] (Ipsi) ; 

if member (' alldep' , [args [3 .. -1] ] ) then 

tt : =t , op(xO) jOpdpsi) ; xx:=tt; pp:=tt; 
elif member ( ' mindep ' , [args [3 . . - 1] ] ) then 

tt:=t; xx:=op(xO); pp:=op(lpsi) ; 
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else tt:=t; xx: =t , op(xO) ; pp :=t , op(lpsi) ; fi: 

if Sr>l then vX:=Vector( [seq(X[i] (xx) , i=l..Sr)]); 

else vX:=Vector([X(xx)]); fi; 
if Sr>l then vPSI :=Vector [row] ( [seqCPSI [i] (pp) , i=l..Sr)]); 

else vPSI:=Vector[row] ([PSKpp)]) ; fi; 
syseqd:=-C vpsi.( map(<iiff ,vphi,t)*T(tt) 

+Matrix( [seq(map(dif f , vphi ,i) , i=xO)] ) . vX 
+vphi*diff (T(tt) ,t)-map(diff ,vX,t) )+vPSI.vphi, 
convert (-vPSI+(vpsi. vphi) *Vector [row] ( [seq(dif f (T(tt) ,i) ,i=xO)]) 
-vpsi .Matrix ( [seq(map(dif f , vX, i) , i=xO)] ) , 'list ' ) [] , 
convert ( (vpsi .vphi) *Vector [row] ( [seq(dif f (T(tt) , i) , i=lpsi)] ) 
-vpsi .Matrix ( [seq(map(diff ,vX, i) , i=lpsi)] ) , 'list')[]} minus {O}: 
lstGerad: = [T(tt) , convert (vX, 'list ' ) [] , convert (vPSI , 'list') []] ; 
if Hi='nohint' then sol : =pdsolve(syseqd, IstGerad) ; 
else sol : =pdsolve (syseqd, IstGerad, HINT=Hi) ; fi; 

if not member ( 'allconst ', [cirgs [3. . -1] ] ) then sol:=reduzConst(sol) ; fi: 
valGerad:=subs (sol, IstGerad) ; 

sol :=[(lstGerad[i]=valGerad[i] )$i=l . .nops (IstGerad)] ; 
sol : =collect (expand(simplif y (sol) ) , [t , op(xO) ,op(lpsi)] ) ; 
if not member (' showdep' , [args [3 .. -1] ] ) then 

xieta: = [xi,seq(eta[i] ,i=l. .n)] ; 

sol:=subs(map(i->i=op(0,i) .IstGerad) ,sol) ; 
else xieta: = [xi(tt) , seq(eta[i] (xx) ,i=l . .n)] ; fi; 
if n=l then xieta:=subs(eta[l]=eta,xieta) fi; 
sol:=subs('PSI'='Psi' , sol); 
sol :=subs(map(i->rhs (i)=lhs (i) ,mapx) , sol) ; 
xieta:=subs(map(i->rhs(i)=lhs(i) ,mapx) ,xieta) ; 
sol2: = [xieta[l] =rhs(sol [1] ) , 

seq(xieta[i+l]=rhs(sol[2+sum(r ['k'] , 'k' = l ..i-1)]), i=l..n)]; 
if member (' split ', [args [3 .. -1] ] ) then sol2 : = [split (sol2)] 
else sol2:=[sol2] fi; 

if member (' showgen' , [args [3. .-!]] ) then sol: = [sol,sol2[]] ; 
else sol:=sol2 fi; 

if n=l then xO : =op(0 ,depvars) else xO:=map(i->op(0,i) .depvars) [] fi; 
sol:=subs({map(i->i(t)=i, [xO,op(lpsi)] ) []}, sol) ; 
if member ( ' showt ' , [args [3 . .-!]]) then 

sol :=subs({map(i->i=i(t) , [xO,op(lpsi)] ) []>,sol) f i; 
return sol [] ; 
end proc: 
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